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ABSTRACT 

We analyzed Chandra observations of three gravitational lenses, 
SBS 0909+523, FBQS 0951+2635, and B 1152+199, to measure the dif- 
ferential X-ray absorption and the dust-to-gas ratio of the lens galaxies. We 
successfully detected the differential X-ray absorption in SBS 0909+523 and 
B 1152+199, and failed to detect it in FBQS 0951+2635 due to the dramatic 
drop in its flux from the ROSAT epoch. These measurements significantly 
increase the sample of dust-to-gas ratio measurements in cosmo logically-distant, 
normal galaxies. Using the larger sample, we obtain an average dust-to-gas 
ratio of E(B — V)/Nn= (1.5 ±0.5)x 10 -22 mag cm 2 atoms -1 with an estimated 
intrinsic dispersion in the ratio of ~40%. This average dust-to-gas ratio is 
consistent with our previous measurement, and the average Galactic value of 
1.7 x 10~ 22 mag cm 2 atoms -1 and the estimated intrinsic dispersion is also 
consistent with the 30% observed in the Galaxy. A larger sample size is still 
needed to improve the measurements and to begin studying the evolution in the 
ratio with cosmic time. We also detected X-ray microlensing in SBS 0909+523 
and significant X-ray variability in FBQS 0951+2635. 

Subject headings: galaxy: evolution - ISM: dust-to-gas ratio 



1. Introduction 



The inter-stellar medium (ISM) and inter-galactic medium (IGM) are ubiquitous, and 
understanding the ISM is crucial to many areas of astronomy such as cosmology (e.g., us- 
ing SNe la or GRBs as cosmological probes), the evolution of galaxies, and star forma- 
tion (e.g., Fall et al. 1996; Madau et al. 1998; Schneider et al. 2004; Corasaniti 2006; 
Ghirlanda et al. 2006; Lacey et al. 2008). The dust-to-gas ratio is a basic property of the 
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ISM, and it is difficult to measure in distant galaxies. In the Galaxy, it is measured by 
comparing the optical extinction and Ly Q absorption towards stars, and a typical value is 
E(B — V)/Nn = 1-7 x 10~ 22 m ag cm 2 atoms -1 with a scatter about the mean of order 30% 



( IBohlin. Savage & Drake Ill978l ). Subsequent studies have extended these measurements to 
additional lines of sight in our Galaxy and to the LMC and SMC, but the basic results are 
little altered aside from some evidence that the dust-to-gas ratios of the LMC and SMC are 
modestly lower (~ 30%) than in the Galaxy (see the review by Draine 2003 and references 
therein). This absorption method is hard to apply to extragalactic sources where we lack 
precise knowledge of the intrinsic source spectra and tend to be probing low density parts of 
the IGM using absorption lines in the spectra of quasars or 7-ray bursts (e.g., Kacprzak et 
al. 2008) rather than the inner regions of normal galaxies. An alternate approach is to use 
emission by the dust in the far-infrared to estimate a dust mass and then compare it to the 
estimated gas mass (e.g., Contini & Contini 2007). 

Gravitational lensing provides a means of studying the properties of the ISM through 
differential absorption in the lens galaxies between the locations of the quasar images. By 
studying the differences in the absorption (or extinction) between multiple images, any 
contamination from the Galaxy in the foreground (Milky Way) or the quasar host galaxy 
in the background is essentially eliminated and we can be confident that we are probing 
the ISM of the lens galaxy. Moreover, since we are comparing two images of the same 
quasar, the method has many of the quantitative advantages of the local approach using 
stars of known spectral type. In optical bands, this technique has been widely used to 
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Ehasdottir et al.ll2006l ). In X-rays, iDai et al. 



1991 



2003 



( 120031 ) and lDai fc Kochanek I (j2005l ) measured the differential c olumn densities in two lenses 
Q 2237+0305 and B 1600+434. Subsequently, bai et all j2006h expanded the results to four 
lenses and proposed that the evolution of the dust-to-gas ratio can be determined given a 
larger sample. In this paper, we present new measurements based on Chandra observations 
of the quasar lenses SBS 0909+523 (Kochanek et al. 1997), FBQS 0951+2635 (Schechter et 
al. 1998), and B 1152+199 (Myers et al. 1999). 



2. Observations and Data Reduction 



We observed SB S 0909+523 and FBQ S 0951+2635 with t he Advanced CCD Ima ging 
Spectrometer (ACIS. iGarmire et al. 1 120031 ) on board Chandra fjWeisskopf et al. 1 120021 ) for 
19.6 and 34.2 ks on December 17, 2006 and March 24, 2007, respectively. B 1152+199 was 
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observed in a separate program (PI: K. Pedersen) with Chandra-ACIS for 16.4 and 8.3 ks 
on February 22, 2005 and February 26, 2005. The back-side illuminated S3 chip was used 
in all the observations. A journal of the Chandra observations is presented in Table [TJ The 
Chandra data were reduced using the CIAO 4.0 software tools provided by the Chandra 
X-ray Center (CXC), following the standard threads on the CXC websiteJll We used the 
most recently reprocessed data products (Reprocessing III) and calibration files (CALDB 
3.4.2). Only events with standard ASCA grades of 0, 2, 3, 4, and 6 were used in the analysis. 
We improved the image quality of the data by removing the pixel randomization applied to 
the event positions by the standard pipeline. In ad dition, we applied a subpixel resolution 



technique (ITsunemi et al. I l200ll ; iMori et al. 1 1200 ll ) to the events on the S3 chip of ACIS 



where the lensed images are located. Figure [T] shows the resulting images. 



3. Spectral Analysis 



We fit the spectra of the lensed quasars using XSPEC VI 1.3.1 ( Arnaudlll996l ) over the 



0.3-8 keV observed energy range. The X-ray spectrum of the ith image was modeled as 

Ni(E,t) = N 0ti (t-Ati)l — \ 

exp {-a(E)N HtGal - a [E(l + *,)] N n>i - a [E{1 + z s )} N H , Src (t - At,)} ,(1) 

where Ni(E,t) is the number of photons per unit energy interval, Nn Gal , Nn ti , and Nn Src 
are the equivalent hydrogen column densities in our Galaxy, the lens galaxy at the position 
of image i, and in the source galaxy, respectively, cr(E) is the photo-electric absorption cross- 
section, and Ati is the time-delay. Although, absorption by the Galaxy, lens and source is 



1 The CXC website is at |http://cxc. harvard.edu/ 



Table 1. Chandra Observations of Gravitational Lenses SBS 0909+523, 
FBQS 0951+2635, and B 1152+199. 



Lenses 


Date 


Exp (sec) 


Count Rate (ct/s, 0.3-8 keV) 


Raw Flux Ratio (A/B) a 


SBS 0909+523 


2006-12-17 


19654 


0.162 + 0.003 


0.34 + 0.02 


FBQS 0951+2635 


2007-03-24 


34243 


0.0074 ± 0.0005 


5.2 ± 1.4 


B 1152+199 


2005-02-22 


16449 


0.189 + 0.003 


3.06 ± 0.14 


B 1152+199 


2005-02-26 


8282 


0.198 + 0.005 


3.83 ± 0.26 



a The raw flux ratio has not been corrected for the differential absorption between the two images. 



-4- 



degenerate when fitting the spectrum of each image, the differential absorption A Ah between 
the lensed images is not and can be determined from a simultaneous fit to the spectra of the 
images. The time-delay at the source complicates the measurement, since if the absorption 
at source is variable (e.g., PKS 1830—211, Dai et al. 2008), it will not be exactly canceled 
when comparing the two spectra. However, since the photo-electric absorption cross-section 
decreases with energy a(E) oc E~ 3 , the absorption at the source will affect the spectrum 
less. Thus, we neglected the time-delay between the images and the absorption at the source 
redshift in this analysis. We note that a significant change in the spectral index on the time 
scale of the time delays will affect the results, especially for quasars with poor S/N spectra 
where the spectral index is degenerate with the absorption. The problem is less severe for 
quasars with moderate to high S/N spectra where the spectral index is less degenerate with 
A^h absorption during spectral modeling. 

A second systematic error we must consider is microlensing of the quasar by the stars 
in the lens galaxy. Our analysis is only affected by chromatic microlensing, where there 
are significant changes in the optical and X-ray spectra that alter the estimated extinction 
or absorption. Indeed, it needs only be achromatic "locally" to separate optical and X-ray 
wavelength/energy regimes, as global shifts of the optical spectrum relative to the X-ray 
spectrum by microlensing also have no effect on the results. For the systems we consider, 
we have only limited control on the level of microlensing. While X-ray microlensing has 
been detected (Chartas et al. 2002, 2007; Dai et al. 2003; Blackburne et al. 2006; Pooley 
et al. 2007), X-ray spectral changes due to microlensing have yet to be detected in large 
part due to the limited sensitivity of the data. Chromatic microlensing in the optical has 
been observed (Poindexter et al. 2007; Anguita et al. 2008; Eigenbrod et al. 2008) and can 
create systematic problems for estimates of the extinction as we had already discovered for 
HE 1104—1805 (Dai et al. 2006). We can, however, compare the optical and X-ray flux 
ratios, since the evidence to date is that the X-ray sources are much more compact than the 
optical (Kochanek et al. 2006; Morgan et al. 2008; Chartas et al. 2008) and so more strongly 
affected by microlensing. If the extinction and absorption-corrected flux ratios of the images 
agree, this is a good indication that our estimates are little affected by microlensing. 

We used the standard wabs and zwabs models in XSPEC to model the Galactic ab- 
sorption a nd absorption at the lens redshi ft. The wabs and zwabs models use cross sec- 



tions from iMorrison fc McCammon I (119831 ) and assume a solar element al abundance from 
Anders fc Ebihara ( 19821 ). The Galactic Ah is fixed using the values from Dickey fc Lockman 



(119901 ) . We fixed the intrinsic power-law index T to be the same for both lensed images, and 
allowed the absorption at the lens to differ. Since the average lens absorption and the power- 
law index are correlated in the spectral modeling, the uncertainties in absorption by the lens 
for the individual images are correlated. To accurately measure the uncertainties in the 
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differential absorption, we used XSPEC to explore the full parameter space of the lens column 
densities and the power-law index. The fitting results are listed in Table [21 and the spectra 
and best fit models are shown in Figures [2J [31 and [U We also obtained the unabsorbed flux 
ratios between the quasar images as part of the model, and verified that adding absorption 
in the source does not alter the estimates of the differential column densities. We describe 
the results for each of the three lenses in the following sub-sections. 



3.1. SBS 0909+523 

SBS 0909+523 (Kochanek et al. 1997; Oscoz et al. 1997) is a two-image lens system 
with a source redshift of z s = 1.377 and a lens redshift of z t = 0.83. SBS 0909+523 is 
well studied in the optical bands, and the extinction curve has been measured using the 
differential extinction method (Mediavilla et al. 2005) to be a Galactic extinction curve with 
a strong 2175A extinction feature (Motta et al. 2002). 

It is clear from the Chandra spectra that image B is more heavily absorbed than image 
A. We estimate that the difference in the column densities is ANu_b,a — (0-0551o;q22) x 
10 22 cm -2 . There are several differential extinction measurements for this system (Falco et al. 
1999; Motta et al. 2002; Mediavilla et al. 2005), and we used the most recent measurement, 
AE(B - V) = 0.32 ± 0.01, from Mediavilla et al. (2005). Combined with our differential 
absorption measurement, we obtained a dust-to-gas ratio of E(B — V)/Nn = (5.8jJ; 7 ) x 10~ 22 
mag cm 2 atoms -1 . 

The optical flux ratio of SBS 0909+523 has changed little from its discovery in 1997 
through 2006 (Kochanek et al. 1997; Motta et al. 2002; Goicoechea et al. 2008) suggesting 
that there is little microlensing variability even if there is evidence for microlensing from the 
difference between the continuum and emission line flux ratios (Motta et al. 2002; Mediavilla 
et al. 2005). The extinction and absorption corrected optical (A/B = 0.35 + 0.02, Mediavilla 
et al. 2005) and X-ray flux ratios {A/B = 0.32 ± 0.03) are consistent but differ from the 
if -band (A/B = 1.12, Lehar et al. 2000) and extinction-corrected emission line flux ratios 
(Mediavilla et al. 2005), indicating that there is microlensing in this system but it is affecting 
the optical and X-ray similarly. The lack of variability and the similar optical/X-ray flux 
ratios suggests the lens is in a region of relatively uniform microlensing magnification. We 
also note that the X-ray variability of the source is modest. SBS 0909+523 was observed three 
times by ROSAT, and the source showed moderate X-ray variability of about 40% (Chartas 
2000). The flux observed by Chandra is in between the three ROSAT flux measurements, 
suggesting no significant total X-ray variability. 



Table 2. Spectral Fitting Results 



Quasar 



Model* 



Galactic Nn N h (B - A) 
(10 22 cm- 2 ) (10 22 cm" 2 ) 



r 



Flux Ratio (A/B) 



X 2 (dof) 



C-Stat(dof) 



FBQS 0951+2635 



SBS 0909+523 



wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 
wabs (zwabs (pow) ) 



0.0169 (fixed) 
0.0229 (fixed) 
0.0273(fixed) 
0.0273(fixed) 
0.0273(fixed) 




0.48 ± 0.06 
0.48 + 0.08 
0.48 ± 0.04 



1.63 + 0.05 
1.32 (fixed) b 
2.08 ± 0.05 
2.10 ± 0.05 
2.07 + 0.03 



0.32 ± 0.03 
3.5 ± 1.6 
1.87 + 0.14 
2.28 + 0.25 



124.5(148) 



56.9(29) 



B 1152+199 (Epoch I) 
B 1152+199 (Epoch II) 
B 1152+199 (Combined) 



148.23(136) 

54.8(75) 
213.1(217) 



Note. — Simultaneous fits to the individual spectra of the lensed quasars. The spectra are constrained to have the same intrinsic power-law 
photon index but can have different absorption column densities at the redshift of the lens. 

a Wabs, zwabs, and pow are the XSPEC models used for the Galactic absorption, absorption at the lens, and the intrinsic power-law spectrum, 
respectively. 

b Thc photon index is fix to the best fit value performed in the 2-8 keV. 
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3.2. FBQS 0951+2635 

FBQS 0951+2635 is a two-image lens discovered by Schechter et al. (1998) with a 
source redshift of z s = 1.24 and a lens redshift of zi = 0.26 (Eigenbrod et al. 2007). 
FBQS 0951+2635 was detected in the ROSAT All Sky Survey (RASS) with a count rate 
of 0.02 ct s _1 in the 0.1-2.4 keV band (Voges et al. 2000). Our Chandra image shows 
that FBQS 0951+2635 is the brightest source within a 2' radius circle, which indicates that 
the source detected in RASS is FBQS 0951+2635. However, the Chandra count rate of 
FBQS 0951+2635, 0.007 ct s _1 in the 0.3-8 keV band, is significantly lower than the pre- 
dicted count rate (0.06 ct s -1 ) from the RASS, suggesting significant source variability. 

The unexpected drop in flux resulted in low S/N spectra that cause significant challenges 
in measuring the differential absorption. We used the C-statistic (Cash 1979), suitable for 
small number statistics, to fit the spectra. Using our standard technique, we found a power- 
law index of F = 0.81 ±0.15 for the source. We failed to detect the differential absorption in 
this lens, as the differential absorption between the two images is consistent with zero. Note, 
however, that the upper limit for the lens absorption in image B, Nn B < 0.22 x 10 22 cm -2 
is larger than that for image A, Aha < 0.05 x 10 22 cm -2 , and that the best fit power-law 
index (r = 0.81 ± 0.15) is unusually flat (e.g., Reeves & Turner 2000; Saez et al. 2008). 
A more plausible interpretation is that the low energy part of the spectrum is absorbed. 
Since the absorption is correlated with the power-law index in the spectral fitting, which is 
a severe problem for low S/N spectra, as an experiment we first fit the spectra in the hard 
band (2-8 keV) to find a steeper, and more sensible, intrinsic power-law (r = 1.32). Then 
we held this fixed and fit for the absorption to find a differential absorption of ANn B A = 
(0.49^o;4i) x 10 22 cm -2 . Due to the low significance of this detection (or non-detection at 
all) because of the low S/N spectra, we will be unable to include this lens in an estimate of 
the mean dust-to-gas ratio. 



3.3. B 1152+199 

B 1152+199 (Myers et al. 1999) consists of two z s = 1.019 quasar images lensed 
by a zi = 0.439 galaxy. The dust extinction curve of B 1152+199 was measured to be 
1.3 < R v < 2.0 (Toft et al. 2000) and R v = 2.1 + 0.1 (Ehasdottir et al. 2006). The 
Ehasdottir et al. (2006) results were based on improved data so we adopt their results. The 
differential absorption can be clearly seen by comparing the Chandra spectra of the two 
images (Figure H]). Since B 1152+199 was observed by Chandra twice, we used our standard 
method and simultaneously fit the spectra for both epochs. We detected a large differen- 
tial absorption AA H b,a = (0.48 ± 0.04) x 10 22 cm -2 . We also fit the two Chandra epochs 
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Fig. 1.— Chandra images of SBS 0909+523, FBQS 0951+2635, and B 1152+199. 
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Fig. 2. — Spectra of images A (bottom) and B (top) of SBS 0909+523. The lower panel 
shows the contributions of the residuals to the \ 2 fit statistics. 
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Fig. 3. — Spectra of images A (top) and B (bottom) of FBQS 0951+2635. 
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Fig. 4. — Spectra of images A (top) and B (bottom) of B 1152+199 for both Chandra 
observations. The black and red lines are the fits for the spectra A and B in the first 
Chandra epoch, and the green and blue lines are for the second epoch. The lower panel 
shows the contributions of the residuals to the x 2 fit statistics. 
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separately, and obtained consistent differential absorptions for both observations (Table [2]). 
We used the model-averaged extinction of E(B — V) — 1.20 ± 0.05 mag from Ehasdottir 
et al. (2006), where the small scatter in the value between their models would have little 
effect on our results. Combining the differential extinction and absorption measurements, 
we found a dust-to-gas ratio for the lens galaxy of E(B — V)/Nu = (2.5 ± 0.2) x 10 -22 mag 
cm 2 atoms -1 . This is the most accurate dust-to-gas ratio measurement we have obtained so 
far. There is weak evidence (l-4er) for a change in the X-ray flux ratios, where the flux ratio 
(A/B) is larger in the second epoch (see Table [2]) • The X-ray flux ratios after correcting for 
absorption (A/B = 1.87±0.14 and 2.28 ±0.25) are also consistent with the range of optical 
flux ratios after correcting for extinction (e.g., AM = 0.85 ± 0.07, 0.6 ± 0.1, Ehasdottir et 
al. 2006). 



Discussion 



We summarize all available differential absorption measurements in Table El including 
the three results from this paper (SBS 0909+523, FBQS 0951+2635, and B 1152+199) and 
the five existing measurements from Dai et al. (2006). We find an average dust-to-gas ratio 
of 

E(B - V)/N H = (1.5 ± 0.5) x 10 -22 mag cm 2 atoms -1 (2) 

with an intrinsic scatter of ~ 40%, where we must exclude FBQS 0951+2635 and HE 1104—1805 
from the fits. If we allow for no intrinsic scatter, we found a ratio of (2.3 ± 0.2) x 10 -22 
mag cm 2 atoms -1 with x 2 /N do f = 3.0, indicating the existence of the additional scatter we 
include in our standard estimate. These estimates are consistent with our previous esti- 
mate in Dai et al. (2006) and t he average Galactic value of 1.7 x 10 -22 mag cm 2 atoms -1 
( IBohlin. Savage fc Drake Ill978l ). The estimates for the intrinsic scatter are also consistent 
with local estimates given our sm all sample size. We compare the results for the lens galaxies 
to the local stellar sample from IBohlin. Savage fc Drake I (119781 ) in Figure [5j and plot the 
evolution of dust-to-gas ratios in Figure [6J Contini & Contini (2007) examined this question 
by modeling the bulk emission of luminous IR galaxies, finding that the dust-to-gas ratio is 
significantly higher than the Galactic value in star forming and central nucleus regions (see 
also Maiolino et al. 2001a,b) and consistent with the Galactic value elsewhere. However, 
the scatter of their results is significantly larger than our more direct measurements. In 
addition, Menard & Chelouche (2008) found several Mgll absorbers at z ~ 1, which also 
have dust-to-gas ratios close to the Galactic value. In our present results, we see no evidence 
for evolution in the dust-to-gas ratio to z ~ 1, an d with a larger sample we will be able t o 
test the ISM evolution models of simulations (e.g.. lDwek II 19981 : Edmunds! l200ll ; Ilnoueil2003l ). 
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Fig. 5. — ATVh versus AE(B — V) for the gravitational lenses (squares and crosses) and Nn 
versus E(B — V) for Galactic stars from Bohlin et al. (1978, filled circles). The dashed line 
shows the mean Galactic relation. The two crosses are the new Chandra measurements for 
SBS 0909+523 and B 1152+199. 
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Fig. 6. — The the dust-to-gas ratio (squares) versus redshift in Galactic units of 1.7 x 1CT 22 
mag cm 2 atoms -1 . The dashed and dotted lines show the Galactic value and its intrinsic 
scatter of 30%. 
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The primary systematic uncertainty for our results comes from chromatic microlensing 
of the quasar by the stars in the lens galaxy distorting the optical and X-ray spectra. Given 
the available information for the systems we consider, it is possible that HE 1104—1805 and 
FBQS 0951+2635 are affected by this problem while the rest show no evidence of a prob- 
lem. In addition, since we have confined our studies to systems with significant evidence 
for extinction, we are relatively safe from minor chromatic microlensing. If microlensing is 
creating the large color differences between the images then we should also see rapid changes 
in the colors because a strong chromatic effect requires the source to lie near or on a caustic. 
Nonetheless, it is possible from two of the systems (HE 1104-1805 and FBQS 0951+2635) 
that microlensing can mimic extinctions of order AE(B — V) = 0.1 mag. However, that the 
lenses follow the local dust-to-gas ratio and with similar scatter is strong evidence that the 
remaining systems are little affected by high levels of chromatic microlensing effect. Monitor- 
ing these systems at both optical and X-ray wavelengths would allow secure determination of 
the absorption properties while simultaneously allowing us to determine the spatial structure 
of the source quasar. 

We thank A. Eliasdottir, B. Menard, L. J. Goicoechea, and the anonymous referee for 
helpful comments. We gratefully acknowledge the financial support by CXC grant G07-8080. 
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Table 3. The Dust-To-Gas Ratio of High Redshift (z > 0) Galaxies 



AJV H AE(B - V) E(B - V)/N H 

Lens z l type Between Images (10 22 cm -2 ) (mag) (10~ 22 mag cm 2 atoms -1 ) 



New Measurements 



SBS 0909+523 0.83 elliptical B, A ^-^-0022 °- 32 ± °- 01 5 - 8 -3 7 

FBQS 0951+2635 0.26 elliptical B, A OA^l'.it -0.12 ± 0.02 • ■ • a 

B 1152+199 0.439 elliptical B, A 0.48 + 0.04 1.20 + 0.05 2.5 + 0.2 



Dai ct al. (2006) Sample 



MG 0414+0534 


0.9584 


elliptical 


A, B 


0.33 + 0.10 


0.18 + 0.11 


0.55 + 0.33 


HE 1104-1805 


0.729 


elliptical 


A, B 


0.055 ± 0.030 


-0.07 + 0.01 


a 


B 1600+434 


0.41 


spiral 


B, A 


o.26i8:S 


0.10 ± 0.03 


0.38 + 0.25 


PKS 1830-211 


0.886 


spiral 


B, A 




3.00 + 0.13 


1.7 + 0.6 


Q 2237+0305 


0.0395 


spiral b 


A, C 


0.04 ± 0.03 


0.11 ± 0.03 


2.8 + 2.2 



a The dust-to-gas ratio was not estimated for FBQS 0951+2635 and HE 1104—1805 because the differential extinction and 
Nn column measurements give opposite signs. 

b Q 2237+0305 is lensed by the central bulge of a spiral galaxy. 
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